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Abstract. Millisecond pulsars with strong magnetic 
fields may be formed through several processes, e.g. 
accretion-induced collapse of magnetized white dwarfs, 
merger of two neutron stars, and accretion-induced phase 
transitions of neutron stars. During the birth of such a 
pulsar, an initial fireball available for a gamma-ray burst 
(GRB) may occur. We here study evolution of a post- 
burst relativistic fireball with energy injection from the 
pulsar through magnetic dipole radiation, and find that 
the magnitude of the optical afterglow from this fireball 
first decreases with time, subsequently flattens, and fi- 
nally declines again. This may provide a natural expla- 
nation for the behavior of the lightcurve of the afterglow 
of GRB970228 if this burst resulted from the birth of a 
strongly magnetic millisecond pulsar. 
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1. Introduction 

The observational results from the BATSE detector on 
CGRO (Fishman & Meegan 1995) strongly suggested that 
the sources of gamma-ray bursts (GRBs) are at cosmolog- 
ical distances or in the galactic halo, leading to a great 
debate on the origin of GRBs. The recent discovery of 
fading sources at X-ray (Costa, et al. 1997a, b; Yoshida, 
et al. 1997) and optical wavebands (Groot, et al. 1997; 
van Paradijs, et al. 1997; Sahu, et al. 1997a, b, c), and the 
detection of the redshift of the counterpart to GRB970508 
(Metzger, et al. 1997), clearly demonstrate that GRBs arc 
at cosmological distances. 
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Most of the models proposed to explain the occur- 
rence of cosmological GRBs are related with compact ob- 
jects from which a large amount of energy available for 
an extremely relativistic fireball is extracted through neu- 
trino emission and/or dissipation of electromagnetic en- 
ergy. One such possibility is newborn millisecond pulsars 
with strong magnetic fields, which may be formed by sev- 
eral models, e.g. (1) accretion-induced collapse of magne- 
tized white dwarfs, (2) merger of two neutron stars, and 
(3) accretion-induced phase transitions of neutron stars. 
In the first model, an accreting white dwarf, when its mass 
increases up to the Chandrasekhar limit, may collapse to 
a rapidly rotating neutron star, and the magnetic field of 
the neutron star may become very strong due to magnetic- 
flux conservation (Usov 1992) or efficient dynamo action 
(Duncan & Thompson 1992). During the birth of the neu- 
tron star, an initial fireball may occur through neutrino 
processes (Dar, et al. 1992) and/or electromagnetic pro- 
cesses (Usov 1992; Yi & Blackman 1997). 

Second, it is usually thought that a neutron-star bi- 
nary will eventually merge into a black hole due to gravi- 
tational radiation. However, since the equations of state at 
high densities are likely stiff for several observational and 
theoretical facts summarized by Cheng & Dai (1997), Dai 
& Cheng (1997), and Cheng & Dai (1998), neutron-star 
binaries like the Hulse- Taylor system might merge into 
massive neutron stars which have both millisecond peri- 
ods and very strong magnetic fields generated by some 
physical processes, e.g. dynamo action powered by tidal 
heating (Vietri 1996). During the coalescence and merging 
of two neutron stars, an initial fireball may be produced 
through neutrino-antineutrino annihilation (Eichler, et al. 
1989; Mathews, et al. 1997) and/or efficient energy trans- 
fer from the orbital energy (Vietri 1996). Third, it has 
been proposed by Cheng & Dai (1996) that some neutron 
stars in low-mass X-ray binaries accrete sufficient mass to 
undergo phase transitions to become strange stars, and af- 
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ter the phase transitions fireballs of GRBs naturally occur 
from the quark surfaces of the hot strange stars (Cheng 
& Dai 1996; Usov 1998) because of very-low-mass crusts. 
The resulting strange stars may have short periods due 
to accretion-induced spin-up prior to the phase transi- 
tions, and they could have strong magnetic fields (Cheng 
& Dai 1997). Therefore, it can be seen that two common 
products of these three energy-source models are newborn 
millisecond pulsars with strong magnetic fields and initial 
fireballs available for GRBs. Furthermore, Kluzniak & Ru- 
derman (1997) proposed that strongly magnetized neutron 
stars differentially rotating at millisecond periods may be 
the central engine of gamma-ray bursts. 

It is natural to expect that if a GRB results from the 
birth of a strongly magnetic millisecond pulsar, then after 
the main GRB the pulsar continuously supplies energy to 
a relativistic fireball through magnetic dipole radiation. In 
this Letter we study evolution of the radiation from such 
a postburst fireball. 

2. The model 

We assume, as a pulsar is born by the models summarized 
in the introduction, an amount of energy E comparable to 
that observed in gamma rays, E ~ 10 51 ergs, is released 
over a time interval less than 100 seconds. This initial 
fireball will expand and accelerate to relativistic velocity 
because of the huge optical depth in the source (Paczynski 
1986; Goodman 1986; Shcmi & Piran 1990). Perhaps, in- 
ternal shocks are formed due to collisions between the 
shells with different Lorentz factors in this period (Rees 
& Meszaros 1994; Paczynski & Xu 1994). During the ac- 
celeration, the initial radiation energy will be converted 
to bulk kinetic energy of the outflow. Subsequently, the 
expansion of the fireball starts to be significantly influ- 
enced by the swept-up interstellar medium (ISM) and two 
shocks will be formed: a forward blastwave and a reverse 
shock (Rees & Meszaros 1992; Katz 1994). A GRB will 
be produced once the kinetic energy is dissipated and ra- 
diated as gamma rays through synchrotron or possibly 
inverse-Compton emission from the accelerated electrons 
in the shocks (Meszaros, et al. 1994; Sari, et al. 1996). The 
postburst fireball will be decelerated continuously, and an 
X-ray, optical and/or radio afterglow will be formed, as 
predicted early by many authors (Paczynski & Rhoads 
1993; Katz 1994; Victri 1997b; Meszaros & Rees 1997). 

At the center of the fireball, the pulsar loses its ro- 
tational energy through magnetic dipole radiation, whose 
power is given by 

= 4 x lO^ergss-'BluP-X, (1) 

where B±^ 2 = B s sin#/10 12 G, B s is the surface dipole 
field strength, 9 is the angle between the magnetic dipole 



moment and rotation axis, P ms is the period in units of 1 
ms, and Rq is the stellar radius in units of 10 6 cm. This 
power varies with time as 

L{t)^{l+t/T)-\ (2) 

where t is one measure of time in the burster's rest frame, 
and T is the initial spin-down timescale defined by 

T = q = 5 x 10 8 s B ; 2 12 P,L/ 45j R 6 - 6 , (3) 

where P is the rate of period increase due to magnetic 
dipole radiation and J45 is the stellar moment of inertia 
111 units of 10 45 gem 2 . (Note that the stellar spin-down 
timescale in the burster's rest frame is equal to that in 
the observer's frame.) For t < T, L can be thought as a 
constant; but for t » T, L decays as oc t~ 2 . 

The power of the pulsar is radiated away mainly 
through electromagnetic waves with frequency of uo = 
2-k/P. Once the electromagnetic waves propagate in the 
shocked ISM, they will be absorbed if u p > u>, where 
uj p is the plasma frequency of the shocked ISM. Since 
the electron number density of the shocked ISM in the 
burster's rest frame is n e = 4j 2 n (Blandford & McKee 
1976), where 7 is the Lorentz factor of the shocked ISM 
and n is the electron number density of the unshocked 
ISM, for cj ~ 10 4 s _1 the inequality lo p > u is always valid 
if n > 0.0l7~ 2 cm~ 3 . This shows that the electromag- 
netic waves cannot propagate through the shocked ISM 
and by this means energy can be continuously pumped 
from the pulsar into the shocked ISM. This idea is sim- 
ilar to that of Pacini (1967), who first proposed that a 
pulsar can continuously supply energy to its surrounding 
supernova remnant through magnetic dipole radiation. 

In the following we assume that the expansion of the 
postburst blastwave in uniform ISM is relativistic and adi- 
abatic. At a time t, the shocked ISM energy is given by 
i? s h = 47rr 2 (r/47 2 )7 2 e' (Waxman 1997a), where r « ct is 
the blastwave radius and e' = A^ 2 nm p c 2 is the shocked 
ISM energy density in the comoving frame (Blandford & 
McKee 1976). This energy should be equal to the sum of 
E/2 and the energy which the fireball has obtained from 
the pulsar based on energy conservation: 

E 

4:TTnm p c 2 ~f 2 r 3 = — + / (1 - (3)L(t - r/c)dt , (4) 
2 Jo 

where /? = (1 — 1/7 2 ) 1 / 2 and the factor 1/2 arises from the 
fact that about half of the total energy of the initial fire- 
ball has been radiated away during the GRB phase (Sari 
& Piran 1995). The term (1 — /?) accounts for the fact that 
the fireball moves away from the incoming photons, and 
the power L(t — r/c) shows the fact that the radiation ab- 
sorbed by the fireball at time t was emitted at the retarded 
time t — r/c. According to Eq. (2), this power should de- 
crease as L(t-r/c) oc [1 + (t - r/c)/T}- 2 = (l + t^/T)- 2 , 
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where i© is the observed time. Because the Lorentz fac- 
tor of the fireball at the initial evolution stage decays as 
7 oc t~ 3 / 2 , the timescale in which the fireball has ob- 
tained energy ~ E/2 from the pulsar can be estimated 
by A^E/L, which is measured in the burster's rest frame. 
The corresponding observer-frame timescale (r) is equal 
to this timescale divided by 2j 2 (Waxman 1997b), viz., 
r = 2E/L. We require r ~ 6 days and T ^> r for the 
optical afterglow of GRB970228. 

We now analyze evolution of the afterglow from a 
postburst fireball with energy injection from the pulsar 
through magnetic dipole radiation. First, at the initial 
stage of the afterglow, viz., t® -C r, the first term on 
the right hand of Eq. (4) is much larger than the sec- 
ond term. In this case, one expects that the fireball is not 
significantly influenced by the stellar radiation, and thus 
its expansion evolves based on 7 = 324£ , g{ 8 n 1 1//8 i 3 ^ 8 , 
where E51 = i?/10 51 ergs, m — n/lcm~ 3 , and i© is in 
units of 1 second. According to Dai & Lu (1997), therefore, 
the synchrotron flux (X-ray) integrated between 2 and 10 
keV decays as Fx oc £® ax , where ax = 3/2 — 3(3 — p)/16, 
and the synchrotron flux density at some optical band de- 
clines as Sopt oc i® a ° pt , where a op t — 3(p— l)/4. Here p is 
the index of the power-law distribution of the accelerated 
electrons in the shocked ISM. The observations of the af- 
terglow of GRB970228 have given a x = 1.4 ± 0.2 (Costa, 
et al. 1997a, b; Yoshida, et al. 1997) and a opt = 2.1±^| in 
R band (Galama, et al. 1997). In order to fit these values, 
we require p ~ 3. 

Second, for T > i® » r, according to Eq. (4), the 
energy which the fireball has obtained from the pulsar 
is much larger than E/2, and the expansion of the fire- 
ball is significantly affected by the pulsar's radiation. In 
the case with steady energy supply, the Lorentz factor 
of the fireball decays as 7 oc r -1 / 2 oc i® 1 ^ 4 (Blandford 
& McKee 1976). This scaling law can also be found by 
neglecting the term E/2 in Eq. (4). Here we have as- 
sumed that the power radiated from the pulsar doesn't 
vary with time. The comoving-frame equipartition mag- 
netic field decays as B' oc 7 ex i® 1 ^ 4 , and the synchrotron 
break frequency drops in time as v m oc 7 3 £>' oc . At the 
same time, since the comoving electron number density is 
n^, oc 7 oc i® 1 ^ 4 and the comoving width of the emission 
region Ar' ~ r/7 a t^ 4 , according to Meszaros & Rees 
(1997a) and Wijers et al. (1997), the comoving intensity 
I' v oc n' e B' far' oc i^ 4 . So the observed peak flux density 
increases in time based on S Vm oc t%1 5 I' Vm oc i®. Thus, for 
v v m , the observed optical flux density as a function of 
observed time is 

S v = S Vm {v/v m )-WI* <xt%-*' 2 . (5) 

Therefore, the requirement of p ~ 3 inferred from the 
early-time behavior of the afterglow leads to the result 
that the observed R-band flux at the subsequent stage 
may keep relatively constant. 
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Third, for i© ^> T, the power of the pulsar due to mag- 
netic dipole radiation rapidly decreases as L oc i® 2 , and 
thus the fireball is hardly influenced by the stellar radia- 
tion and the optical flux of the afterglow will significantly 
decline with time again. 

3. Constraints on stellar parameters 

The analysis of the lightcurvc in the last section can 
be directly applied to discussion on the afterglow of 
GRB970228. As pointed out by Galama et al. (1997), the 
initial decrease of the optical afterglow can be easily ex- 
plained by the popular fireball model (Meszaros & Rees 
1997; Wijers, et al. 1997; Waxman 1997a; Victri 1997a; 
Tavani 1997; Reichart 1997; Dai & Lu 1997), but the 
subsequent flattening is not consistent with this simple 
model. Further considering the result observed by HST 
on 4 September (Fruchter, et al. 1997), therefore, we can 
see that the magnitude of the R-band brightness of the af- 
terglow of GRB970228 first rapidly decreased with time, 
subsequently flattened, and finally declined again. The 
lightcurve analyzed in the last section is well consistent 
with these observational results. 

If the afterglow of GRB970228 was radiated from 
a relativistic fireball with energy injection from a pul- 
sar through magnetic dipole radiation, we now give con- 
straints on some parameters of this pulsar. Inserting Eq. 
(1) into t = 2E/L, we get the timescale at which the 
fireball has obtained energy ~ E/2: 

r«5xl0 7 si?5iBl 2 12 P 4 s ^ 6 - (6) 

This equation can be further changed into the following 
form: 

B±,i2 « 10El( 2 (r/6d)-^P^R e 3 . (7) 

The R-band magnitude of the optical afterglow dur- 
ing one month following 6 March after GRB970228 didn't 
vary with time (Galama, et al. 1997), which means that 
during this period the stellar radiation significantly af- 
fected the evolution of the fireball. This in fact requires 
one condition: T > 36 days. On the other hand, the V- 
band magnitude of the afterglow on 4 September dropped 
to V = 28.0 ±0.25 which corresponds to the R-band mag- 
nitude R « 27.0 (Fruchter, et al. 1997), showing that for 
f® » 1 month the stellar radiation played no role in the 
expansion of the fireball. This requires another condition: 
T <C 188 days. Furthermore, after comparing the observed 
data with the decay law of the optical flux at the early 
stage or at the late stage (5 op t oc i© 3 ^ 2 ), we can infer that 
the flattening phase of the optical flux must have ended 
around 60 days. This means T < 60 days. From these 
conditions and Eqs. (3) and (7), we obtain 

1.0 < P ms £ 5 \ /2 (r/6d)- 1/2 J45 1/2 < 1-2 • (8) 
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The equations of state at high densities are likely stiff and 
thus the moments of inertia of rapidly rotating neutron 
stars with mass > 1.4M are J45 ~ 2 (Datta 1988; Weber 
& Glendenning 1993). Further adopting R 6 ~ 1, E 5 i ~ 1 
and t ~ 6 days, we can find 1.4 < P ms < 1.7 and 20 < 
B± t i2 < 30. Therefore, we suggest that the central engine 
of GRB970228 could be a strongly magnetic millisecond 
pulsar. 

4. Discussion 

We have studied evolution of a postburst fireball with 
energy injection from a millisecond pulsar with a strong 
magnetic field through magnetic dipole radiation if the 
initial fireball of a GRB has been produced during the 
birth of the pulsar. In the case of r < T or P nls < 
3.2 mslll" 2 E 5 ^ 2 , according to Eqs. (3) and (6), we found 
that the magnitude of the optical afterglow from this fire- 
ball first decreases with time, subsequently flattens, and 
finally declines again. We have assumed a newborn mil- 
lisecond pulsar with a strong magnetic field to be an 
origin of GRB970228, but have not invoked any spe- 
cific model. There may be several mechanisms relating 
the birth of strongly magnetic millisecond pulsars with 
GRBs, e.g. accretion-induced collapse of magnetized white 
dwarfs, merger of neutron-star binaries, and accretion- 
induced phase transitions of neutron stars. One expects 
that one of them could be the mechanism of the birth 
of a strongly magnetic millisecond pulsar if this pulsar is 
the central engine of GRB970228. In the case of r > T 
or P ms > 3.2 ms/lffijj 1 ' 2 , the millisecond pulsar as the 
central engine cannot supply energy to the postburst fire- 
ball effectively enough for the flattening phase to occur. It 
is interesting to note that this feature doesn't depend on 
magnetic field strength, which influences r and T in the 
same way. 

We further assume that some cosmological GRBs re- 
sult from the birth of strongly magnetic millisecond pul- 
sars. If the effective dipole magnetic-field strength of a 
newborn pulsar is extremely high (e.g., B± = B s sm8 ~ 
10 15 G), the stellar rotational energy is dissipated due to 
magnetic dipole radiation in observed time T < 10 2 s. In 
other words, the fireball of a GRB can obtain a large 
amount of energy from the pulsar only in such a short 
timcscale. For t§ 3> T, thus, the expansion of the post- 
burst fireball is hardly affected by the magnetic dipole 
radiation of the pulsar. In the superstrong-magnetic-field 
case, the flattening behavior might not be observable in 
the optical counterparts of GRBs. 
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